Method and system of representing a sound field

ABSTRACT

The invention relates to a method of representing a sound field. The inventive method includes a step involving the acquisition of measurement signals (c n ) which are delivered by acquisition means ( 1 ) comprising one or more simple sensors ( 2   n ) that are exposed to said sound field (P). The invention is characterised in that it comprises: a step involving the determination of encoding filters which are representative of at least the structural characteristics of the aforementioned acquisition means ( 1 ); and a step whereby the measurement signals (c n ) are processed by applying said encoding filters to the signals (c n ), in order to determine a finite number of representative coefficients over time and in the three-dimensional space of the sound field (P), said coefficients being used to produce a representation of the sound field (P) which is essentially independent of the characteristics of the acquisition means ( 1 ).

The present invention relates to a method and a device for representing an acoustic field from signals issued by acquisition means.

Current methods and systems for acquiring and representing sound environments use models based on acquisition means that are physically impracticable, in particular as far as the electro-acoustic and/or structural characteristics of these acquisition means are concerned.

The acquisition means comprise, for example, a set of measuring elements or elementary sensors arranged in specific spatial locations and having intrinsic electro-acoustic acquisition characteristics.

The current systems are limited by the structural characteristics of the acquisition means, such as the physical arrangement and electro-acoustic characteristics of the elementary sensors, and issue degraded representations of the sound environment to be acquired.

The systems subsumed under the term “Ambisonic”, for example, only consider the directions of the source of sounds relative to the centre of the acquisition means comprising a plurality of elementary sensors, which results in the acquisition means being equivalent to a point microphone.

However, the impossibility of positioning all of the elementary sensors at a single point limits the efficiency of these systems.

Furthermore, these systems represent the sound environment by modelling virtual sources, the angular distribution of which around the centre theoretically allows a sound environment of this type to be obtained.

However, the unavailability of elementary sensors having high directivity characteristics limits these systems to a level of representation precision that is commonly known as “order one”, on a mathematical basis known as the basis of spherical harmonics.

In other systems, such as that employing the method and the acquisition device disclosed in patent application No. WO-01-58209, the acquisition is based on the measurement, in a plane, of information that is representative of the sound environment to be acquired.

However, these systems use models based on optimal elementary sensors that are necessarily arranged on a circle and cause significant amplification of the background noise of the sensors.

These systems therefore require sensors of which the intrinsic background noise is extremely low, and are thus impracticable.

Furthermore, in these systems, the sound environment is only described by a bi-dimensional model, which entails a significant and reductive approximation of the real sound characteristics.

It would therefore seem that the representations of sound environments made by the current systems are incomplete and degraded, and that there is no system that allows a faithful representation to be obtained.

The object of the invention is to solve this problem by providing a method and a device issuing a representation of the acoustic field that is substantially independent of the characteristics of the acquisition means.

The present invention relates to a method for representing an acoustic field comprising a step involving the acquisition of measurement signals issued by acquisition means comprising one or more elementary sensors that are exposed to said acoustic field, characterised in that it comprises:

-   -   a step involving the determination of encoding filters that are         representative of at least the structural characteristics of         said acquisition means; and     -   a step involving the processing of said measurement signals by         applying said encoding filters to these signals in order to         determine a finite number of coefficients representative over         time and in the three-dimensional space of said acoustic field,         said coefficients allowing a representation of said acoustic         field to be obtained that is substantially independent of the         characteristics of said acquisition means.

According to other characteristics:

-   -   said structural characteristics comprise at least position         characteristics of said elementary sensors relative to a         predetermined reference point of said acquisition means;     -   encoding filters are also representative of electro-acoustic         characteristics of said acquisition means;     -   said electro-acoustic characteristics comprise at least         characteristics related to the intrinsic electro-acoustic         acquisition capacities of said elementary sensors;     -   coefficients allowing a representation of the acoustic field to         be obtained are what are known as Fourier-Bessel coefficients         and/or linear combinations of Fourier-Bessel coefficients;     -   step involving the determination of the encoding filters         comprises:         -   a sub-step involving the determination of a sampling matrix             that is representative of the acquisition capacities of said             acquisition means;         -   a sub-step involving the determination of an             intercorrelation matrix that is representative of the             similarity between said measurement signals issued by the             elementary sensors forming said acquisition means; and         -   a sub-step involving the determination of an encoding matrix             from said sampling matrix, said intercorrelation matrix and             a parameter that is representative of a desired compromise             between faithfulness of representation of the acoustic field             and minimisation of the background noise caused by the             acquisition means, which matrix is representative of said             encoding filters;     -   sub-steps involving the determination of the matrices are         carried out for a finite number of operating frequencies;     -   step involving the determination of the sampling matrix is         carried out, for each of said elementary sensors forming said         acquisition means, from:         -   parameters that are representative of the position of said             sensor relative to the centre of said acquisition means;             and/or         -   a finite number of coefficients that are representative of             the acquisition capacities of said sensor;     -   step involving the determination of the sampling matrix (B) is         also carried out from at least one of the following parameters:         -   parameters that are representative of the frequency             responses of all or some of the sensors;         -   parameters that are representative of the directivity             patterns of all or some of the sensors;         -   parameters that are representative of the orientations of             all or some of the sensors, i.e. of their maximum             sensitivity direction;         -   parameters that are representative of the power spectral             densities of the background noise of all or some of the             sensors;         -   a parameter specifying the order in which the representation             is conducted;         -   a parameter that is representative of a list of             coefficients, the power of which must be equal to the power             of the corresponding coefficient in the acoustic field to be             represented;     -   it comprises a calibration step allowing all or some of the         parameters used in said step involving the determination of the         encoding filters, to be issued;     -   calibration step comprises, for at least one of said elementary         sensors forming said acquisition means:         -   a sub-step involving the acquisition of signals that are             representative of the acquisition capacities of said at             least one sensor; and         -   a sub-step involving the determination of parameters             representative of electro-acoustic and/or structural             characteristics of said at least one sensor;     -   calibration step further comprises:         -   a sub-step involving the emission of a specific acoustic             field toward said at least one sensor, said acquisition             sub-step corresponding to the acquisition of the signals             issued by this sensor when it is exposed to said specific             acoustic field; and         -   a sub-step involving the modelling of said specific acoustic             field in a finite number of coefficients, in order to allow             said sub-step involving the determination of parameters that             are representative of electro-acoustic and/or structural             characteristics of the sensor to be carried out;     -   said calibration step comprises a sub-step involving the         reception of a finite number of signals that are representative         of the electro-acoustic and structural characteristics of said         sensors forming said acquisition means, which signals are used         directly during said sub-step involving the determination of the         electro-acoustic and/or structural characteristics of said         acquisition means; and     -   it comprises an input step allowing all or some of the         parameters used during said step involving the determination of         the encoding filters, to be determined.

The invention also relates to a computer programme comprising programme code instructions for implementing the steps of the method as described above, when said programme is executed on a computer.

The invention also relates to a movable support of the type comprising at least one operation processor and a non-volatile memory element, characterised in that said memory comprises a programme comprising code instructions for implementing the steps of the method as described above, when said processor executes said programme.

The invention also relates to a device for representing an acoustic field that is connectable to acquisition means comprising one or more elementary sensors issuing measurement signals when they are exposed to said acoustic field, characterised in that it comprises a module for processing the measurement signals by applying encoding filters that are representative of at least the structural characteristics of said acquisition means to these measurement signals, in order to issue a signal that comprises a finite number of coefficients representative over time and in the three-dimensional space of said acoustic field, said coefficients allowing a representation of said acoustic field to be obtained that is substantially independent of the characteristics of said acquisition means.

According to other characteristics of the invention:

-   -   encoding filters are also representative of electro-acoustic         characteristics of said acquisition means;     -   it further comprises means for determining said encoding filters         that are representative of structural and/or electro-acoustic         characteristics of said acquisition means;     -   said means for determining encoding filters receive at the input         at least one of the following parameters:         -   parameters that are representative of the positions,             relative to centre of said acquisition means, of all or some             of the sensors;         -   a finite number of coefficients that are representative of             the acquisition capacities of all or some of the sensors;         -   parameters that are representative of the frequency             responses of all or some of the sensors;         -   parameters that are representative of the directivity             patterns of all or some of the sensors;         -   parameters that are representative of the orientations of             all or some of the sensors, i.e. of their maximum             sensitivity direction;         -   parameters that are representative of the power spectral             densities of the background noise of all or some of the             sensors;         -   a parameter that is representative of the desired compromise             between faithfulness of representation of the acoustic field             and minimisation of the background noise caused by the             acquisition means;         -   a parameter specifying the order in which the encoding is             conducted;         -   a parameter that is representative of a list of             coefficients, the power of which must be equal to the power             of the corresponding coefficient in the acoustic field to be             represented;     -   it is associated with means for determining all or some of the         parameters received by said means for determining the encoding         filters, said means comprising at least one of the following         elements:         -   means for inputting parameters; and/or         -   calibration means;     -   it is associated with means for formatting said measurement         signals, in order to issue a corresponding formatted signal.

A better understanding of the invention will be facilitated by reading the following description, given solely by way of example and with reference to the accompanying drawings, in which:

FIG. 1 is an illustration of a spherical reference figure;

FIG. 2 is a diagram illustrating the acquisition means used;

FIG. 3 is a general flow chart of the method of the invention;

FIG. 4 is a detailed flow chart of an embodiment of the calibration step of the method of the invention;

FIG. 5 is a detailed flow chart of an embodiment of the step involving the determination of the encoding filters of the method of the invention;

FIG. 6 is a detailed diagram of an embodiment of the step involving the application of the encoding filters; and

FIG. 7 is a block diagram of a device that is suitable for carrying out the method of the invention.

FIG. 1 illustrates a conventional spherical reference figure, so as to clarify the coordinate system referred to in the text.

This reference figure is an orthonormal reference figure, having an origin 0 and comprising three axes (OX), (OY) and (OZ).

In this reference figure, a position marked {right arrow over (x)} is described by means of its spherical coordinates (r, θ, φ), wherein r denotes the distance relative to the origin O, θ the orientation in the vertical plane and φ the orientation in the horizontal plane.

In a reference figure of this type, an acoustic field is known if the sound pressure marked p(r, θ, φ, t), the Fourier transform of which is marked P(r, θ, φ, f), wherein f denotes the frequency, is defined at each point and at each instant t.

The method of the invention is based on the use of spatio-temporal functions allowing any acoustic field over time and in three-dimensional space to be described.

In the described embodiments these functions are what are known as spherical Fourier-Bessel functions of the first kind referred to hereinafter as Fourier-Bessel functions.

In a zone devoid of sources and obstacles, the Fourier-Bessel functions correspond to solutions to the wave equation and form a basis that generates all of the acoustic fields produced by sources located outside this zone.

Any three-dimensional acoustic field may thus be expressed by a linear combination of Fourier-Bessel functions, according to the expression of the inverse Fourier-Bessel transform, which is expressed as follows: ${P\left( {r,\theta,\phi,f} \right)} = {4\pi{\sum\limits_{l = 0}^{\infty}{\sum\limits_{m = {- l}}^{l}{{P_{l,m}(f)}j^{l}{j_{l}({kr})}{y_{l}^{m}\left( {\theta,\phi} \right)}}}}}$

In this equation, the terms P_(l,m)(f) are defined as the Fourier-Bessel coefficients of the field p(r, θ, φ, t), ${k = {k = \frac{2\pi\quad f}{c}}},$ c is the velocity of sound in air (340 ms⁻¹), j_(l) (kr) is the spherical Bessel function of the first kind of order l, defined by ${{j_{l}(x)} = {\sqrt{\frac{\pi}{2x}}{J_{l + {1/2}}(x)}}},$ wherein J_(v)(x) is the Bessel function of the first kind of order v, and y_(l) ^(m)(θ, φ) is the real spherical harmonic of order l and term m, with m ranging from −l to l, defined by: y _(l) ^(m)(θ,φ)=P _(l) ^(|m|)(cos θ)trg _(m)(φ) wherein: ${{trg}_{m}(\phi)} = \left\{ \begin{matrix} {\frac{1}{\sqrt{\pi}}{\cos\left( {m\quad\phi} \right)}} & {{{where}\quad m} > 0} \\ \frac{1}{\sqrt{2\pi}} & {{{where}\quad m} = 0} \\ {\frac{1}{\sqrt{\pi}}{\sin\left( {m\quad\phi} \right)}} & {{{where}\quad m} < 0} \end{matrix} \right.$

In this equation, P_(l) ^(m)(x) are the associated Legendre functions, defined by: ${P_{l}^{m}(x)} = {\sqrt{\frac{{2l} + 1}{2}}\sqrt{\frac{\left( {l - m} \right)!}{\left( {l + m} \right)!}}\left( {1 - x^{2}} \right)^{m/2}\frac{\mathbb{d}^{m}}{\mathbb{d}x^{m}}{P_{l}(x)}}$ wherein P_(l)(x) are Legendre polynomials, defined by: ${P_{l}(x)} = {\frac{1\quad\mathbb{d}^{l}}{2^{l}{l!}\quad{\mathbb{d}x^{l}}}\left( {x^{2} - 1} \right)^{l}}$

The Fourier-Bessel coefficients are also expressed in the temporal domain by the coefficients p_(l,m)(t), corresponding to the inverse temporal Fourier transform of the coefficients P_(l,m)(f).

In other embodiments, the acoustic field is decomposed on a function base, wherein each of the functions is expressed by a potentially infinite linear combination of Fourier-Bessel functions.

FIG. 2 illustrates schematically acquisition means 1 comprising N elementary sensors 2 ₁ to 2 _(N).

These elementary sensors are arranged at specific points in space around a predetermined point 4, designated as the centre of the acquisition means 1.

The position of each elementary sensor may thus be expressed in space, in a spherical reference figure such as that described with reference to FIG. 1, centred on the centre 4 of the acquisition means 1.

When exposed to an acoustic field P each sensor 2 _(n) of the acquisition means 1 issues a measurement signal c_(n), which corresponds to the measurement made by the sensor in the acoustic field P.

The acquisition means 1 thus issue a plurality of signals c₁ to c_(N), which are the signals of the measurement of the acoustic field P by the acquisition means 1.

These measurement signals c₁ to c_(N) issued by the acquisition means 1 are thus directly related to the acquisition capacities of the elementary sensors 2 ₁ to 2 _(N).

FIG. 3 illustrates a general flow chart of the method of the invention.

The method starts with a step 10 involving the inputting of parameters and a step 20 involving the calibration of the acquisition means, which allow a set of parameters that are representative of the structural and/or electro-acoustic characteristics of the acquisition means 1 to be defined.

Some parameters, in particular parameters that are representative of electro-acoustic characteristics, are frequency-dependent.

The inputting step 10 and the calibration step 20, which will be described in greater detail with reference to FIG. 4, may be carried out simultaneously or in any order.

Equally, the method of the invention may comprise only the inputting step 10.

The inputting step 10 and the calibration step 20 allow all or some of the following parameters to be determined for one or more sensor:

-   -   parameters {right arrow over (x)}_(n) that are representative of         the position of the sensor 2 _(n) relative to the centre 4 of         the acquisition means 1, which are written in spherical         coordinates (r_(n),θ_(n),φ_(n));     -   parameters d_(n)(f) that are representative of the directivity         diagram of the sensor 2 _(n), which may take any values between         0 and 1 and allows the directivity of the sensor 2 _(n) to be         described by a combination of omnidirectional and bi-directional         diagrams:     -   if d_(n)(f)=0, the sensor is omnidirectional     -   if d_(n)(f)=½, the sensor is cardioid     -   if d_(n)(f)=1, the sensor is bi-directional;     -   parameters α_(n)(f) that are representative of the orientation         of the sensor 2 _(n), i.e. its maximum sensitivity direction,         which is given by the angle couple (θ_(n) ^(α),φ_(n) ^(α))(f);     -   parameters H_(n)(f) that are representative of the frequency         response of the sensor 2 _(n), corresponding, for each frequency         f, to the sensitivity of the sensor 2 _(n) in the direction         α_(n)f);     -   parameters σ² _(n)(f) that are representative of the power         spectral density of the background noise of the sensor 2 _(n);     -   parameters B_(n,l,m)(f) that are representative of the         acquisition capacities of the sensor 2 _(n), i.e. of the manner         in which the sensor 2 _(n) gathers information on the acoustic         field P. Each B_(n,l,m)(f) is thus representative of the         acquisition capacities of a sensor and, in particular, of its         position in space, and the total of B_(n,l,m)(f) is         representative of the sampling of the acoustic field P carried         out by the acquisition means 1;     -   a parameter μ(f) specifying a compromise between faithfulness of         representation of the acoustic field P and minimisation of the         background noise produced by the sensors 2 ₁ to 2 _(N), and         being able to take all values between 0 and 1:         -   if μ(f)=0, the background noise is minimal;         -   if μ(f)=1, the spatial quality is maximal;     -   a parameter L(f) specifying the order in which the         representation is conducted; and     -   a parameter {(l_(k),m_(k))}(f) that is representative of a list         of coefficients, the power of which must be equal to the power         of the corresponding coefficient in the acoustic field to be         represented.

In simplified embodiments, all or some of the described parameters are considered to be frequency-independent.

The parameters μ(f), L(f) and {(l_(k),m_(k))}(f) are representative of optimisation strategies allowing optimal extraction of spatio-temporal information on the acoustic field P from measurement signals c₁ to c_(N), and are inputted during the inputting step 10. The other parameters may be input during the inputting step 10 or determined during the calibration step 20.

In simplified embodiments, the method of the invention is carried out only with the parameters μ(f), L(f) and all of the parameters {right arrow over (x)}_(n), or all of the parameters B_(n,l,m)(f) or a combination of parameters {right arrow over (x)}_(n) and B_(n,l,m)(f), so that there is at least one parameter per elementary sensor 2 _(n).

Of course, all or some of the parameters used may be issued by memories or dedicated devices, it being possible for an operator to equate these processes to the direct inputting step 10, as described.

Following the input step 10 and/or the calibration step 20, the method comprises a step 30 involving the determination of encoding filters that are representative of at least the structural characteristics, and advantageously the electro-acoustic characteristics, of the acquisition means 1.

This step 30, which will be described in greater detail with reference to FIG. 5, allows all the parameters determined during the input step 10 and/or the calibration step 20 to be taken into account.

These encoding filters are therefore representative of at least the position characteristics of the elementary sensors 2 _(n) relative to the reference point 4 of the acquisition means 1.

Advantageously, these filters are also representative of other structural characteristics of the acquisition means 1, such as the orientation or mutual influences of the elementary sensors 2 ₁ to 2 _(N), and also their electro-acoustic acquisition capacities and, in particular, their background noise, their directivity diagram, their frequency response, etc.

The encoding filters obtained at the end of the step 30 may be stored, so that the steps 10, 20 and 30 are only repeated in the event of modification of the acquisition means 1 or optimisation strategies.

These encoding filters are applied during a step 40 involving the processing of signals c₁ to c_(N) derived from the elementary sensors 2 ₁ to 2 _(N).

The processing entails filtering the signals and combining the filtered signals.

Following this step 40 involving the processing of the measurement signals by applying encoding filters thereto, a finite number of coefficients representatives over time and in the three-dimensional space of the acoustic field P is issued.

These coefficients are what are known as Fourier-Bessel coefficients, marked P_(l,m)(f) and correspond to a representation of the acoustic field P that is substantially independent of the characteristics of the acquisition means 1.

It would therefore appear that the method of the invention allows a faithful representation of the acoustic field of which the temporal and spatial characteristics are being transcribed, whatever acquisition means are used.

FIG. 4 illustrates a flow chart of an embodiment of the calibration step 20.

In this embodiment, the calibration step 20 allows the coefficients B_(n,l,m)(f) which are representative of the acquisition capacities of the acquisition means 1, to be determined directly.

This step 20 starts with a sub-step 22 involving the emission of a specific acoustic field toward the acquisition means 1, and with a sub-step 24 involving the acquisition of measurement signals by the acquisition means 1 exposed to the emitted acoustic field.

Theses sub-steps 22 and 24 are repeated for a plurality Q of specific different fields, and require means for generating specific acoustic fields and means for displacing and/or rotating the acquisition means 1.

For example, the calibration step 20 is carried out using means for generating an acoustic field that merely comprise a fixed loudspeaker, which is assumed to be a point loudspeaker having a flat frequency response, the loudspeaker and the acquisition mans 1 being placed in an anechoic environment.

In each generating sub-step 22, the loudspeaker emits the same acoustic field and the acquisition means 1 are placed in the same position, but they are oriented in different and known directions.

It is, of course, also possible to displace the loudspeaker.

Therefore, in the reference figure of the acquisition means 1, the loudspeaker is in a different position (r_(q) ^(hp),θ_(q) ^(hp),φ_(q) ^(hp)) for each field q generated.

The acquisition means 1 are thus exposed to an acoustic field q, the Fourier-Bessel coefficients of which P_(l,m,q)(f), in the reference figure of the acquisition means 1, are known up to a given order, marked L₃.

In the described embodiment, the measurement signals issued following the acquisition sub-step 24 are a finite number of coefficients that are representative of the generated acoustic field q, as well as of the acquisition capacities of the acquisition means 1.

The parameters L₃ and Q are selected so as to respect the condition: Q≧(L₃+1)²

Advantageously, the method subsequently comprises a modelling sub-step 26, allowing a representation of the Q acoustic fields emitted during the sub-step 22 to be determined.

A modelling matrix P that is representative of all of the known fields Q to which the acquisition means 1 are exposed in succession, is thus determined during the sub-step 26. This matrix P is a matrix of the size (L₃+1)² over Q, comprising elements P_(l,m,q)(f), the indices (l,m) designating the row (l²+l+m), and the index q designating the column q. The matrix P therefore has the following form: $\begin{bmatrix} {P_{0,0,1}(f)} & {P_{0,0,2}(f)} & \cdots & {P_{0,0,Q}(f)} \\ {P_{1,{- 1},1}(f)} & {P_{1,{- 1},2}(f)} & \cdots & {P_{1,{- 1},Q}(f)} \\ {P_{1,0,1}(f)} & {P_{1,0,2}(f)} & \cdots & {P_{1,0,Q}(f)} \\ {P_{1,1,1}(f)} & {P_{1,1,2}(f)} & \cdots & {P_{1,1,Q}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {P_{L_{3},{- L_{3}},1}(f)} & {P_{L_{3},{- L_{3}},2}(f)} & \cdots & {P_{L_{3},{- L_{3}},Q}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {P_{L_{3},0,1}(f)} & {P_{L_{3},0,2}(f)} & \cdots & {P_{L_{3},0,Q}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {P_{L_{3},L_{3},1}(f)} & {P_{L_{3},L_{3},2}(f)} & \cdots & {P_{L_{3},L_{3},Q}(f)} \end{bmatrix}$

In the described embodiment, the acoustic field produced by the loudspeaker is modelled by spherical radiation, such that, in the reference figure of the acquisition means 1, the coefficients P_(l,m,q)(f) of each acoustic field q thus generated are known, owing to the relationship: ${P_{l,m,q}(f)} = {\frac{1}{r_{q}^{hp}}{\mathbb{e}}^{- \frac{{j2\pi}\quad r_{q}^{hp}f}{c}}{\xi_{l}\left( {r_{q}^{hp},f} \right)}{y_{l}^{m}\left( {\theta_{q}^{hp},\phi_{q}^{hp}} \right)}}$ wherein ${\xi_{l}\left( {r_{q}^{hp},f} \right)} = {\sum\limits_{k = 0}^{l}{\frac{\left( {l + k} \right)!}{2^{k}{k!}{\left( {l - k} \right)!}}\left( \frac{{j2\pi}\quad r_{q}^{hp}f}{c} \right)^{- k}}}$

The coefficients obtained in the sub-step 26 are then used in a sub-step 28, in order to determine parameters that are representative of structural and/or sound characteristics of the acquisition means 1.

In the described embodiment, this sub-step 28 also uses the modelling matrix P determined in the sub-step 26.

This sub-step 28 starts with the determination of a matrix C that is representative of all of the signals c_(n,q)(t) picked up at the output of N sensors in response to Q known fields. This matrix C is a matrix of size N over Q, comprising elements C_(n,q)(f), the index n designating the row n, and the index q designating the column q. The elements C_(n,q)(f) are deduced from the signals c_(n,q)(t) by Fourier transformation. The matrix C therefore has the following form: $\begin{bmatrix} {C_{1,1}(f)} & {C_{1,1}(f)} & \cdots & {C_{1,1}(f)} \\ {C_{2,1}(f)} & {C_{2,2}(f)} & \cdots & {C_{2,Q}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {C_{N,1}(f)} & {C_{N,2}(f)} & \cdots & {C_{N,Q}(f)} \end{bmatrix}$

The matrix C is representative of the acquisition capacities of the acquisition means 1 and the Q emitted acoustic fields.

In the described embodiment, the coefficients B_(n,l,m)(f) are determined from the matrices C and B, during the sub-step 28, using conventional methods of general matrix inversion, applied to the relationship that links C to P. For example, the coefficients B_(n,l,m)(f) are placed in a matrix B that is determined by the following relationship: B=C P ^(T)(P P ^(T))⁻¹ B is a matrix of size N over (L₃+1)² comprising coefficients B_(n,l,m)(f), the index n designating the row n and the indices (l,m) designating the column l²+l+m. The matrix B therefore has the following form: $\left\lbrack \left. \quad\begin{matrix} {B_{1,0,0}(f)} & {B_{1,1,{- 1}}(f)} & {B_{1,1,0}(f)} & {B_{1,1,1}(f)} & \ldots & {B_{1,L_{3},{- L_{3}}}(f)} & \ldots & {B_{1,L_{3},0}(f)} & \ldots & {B_{1,L_{3},L_{3}}(f)} \\ {B_{2,0,0}(f)} & {B_{2,1,{- 1}}(f)} & {B_{2,1,0}(f)} & {B_{2,1,1}(f)} & \ldots & {B_{2,L_{3},{- L_{3}}}(f)} & \ldots & {B_{2,L_{3},0}(f)} & \ldots & {B_{2,L_{3},L_{3}}(f)} \\ \vdots & \vdots & \vdots & \vdots & \quad & \vdots & \quad & \vdots & \quad & \vdots \\ {B_{N,0,0}(f)} & {B_{N,1,{- 1}}(f)} & {B_{N,1,0}(f)} & {B_{N,1,1}(f)} & \ldots & {B_{N,L_{3},{- L_{3}}}(f)} & \ldots & {B_{N,L_{3},0}(f)} & \ldots & {B_{N,L_{3},L_{3}}(f)} \end{matrix} \right\rbrack \right.$

These sub-steps 26 and 28 are carried out for each operating frequency, and the coefficients thus determined directly form the parameters that are representative of the acquisition capacities of the acquisition means 1.

The sub-steps 26 and 28 of the calibration step 20 may be carried out in various ways, as a function of the parameters that have to be determined.

For example, in the case where the calibration step 20 allows the position {right arrow over (x)}_(N) of each sensor 2 _(n) to be determined, the sub-steps 26 and 28 use the propagation times of the waves emitted by the loudspeakers to reach the sensors 2 _(n). The position of each sensor 2 _(n) is determined using at least three propagation time measurements, according to triangulation methods.

In another case, when the loudspeaker emits a given impulse, the sub-steps 26 and 28 allow the impulse responses of each sensor 2 _(n) to be determined from the signals c_(n,q)(t).

Standard methods for determining impulse responses, such as MLS (maximum length sequence), for example, are used in this case.

Advantageously, the calibration step 20 allows electro-acoustic characteristics of the sensors to be determined. It then starts by determining the directivity diagram of each sensor 2 _(n) for each given frequency f, for example, by determining the frequency response of each sensor 2 _(n) for a plurality of directions.

In a second stage, all or some of the following parameters are determined:

-   -   parameters α_(n)(f) that are representative of the orientation         of each sensor 2 _(n), i.e. of its maximum sensitivity         direction, given by the angles (θ_(n) ^(α),φ_(n) ^(α))(f), for         which the directivity diagram admits a maximum to the common         frequency f;     -   parameters H_(n)(f) that are representative of the frequency         response of each sensor 2 _(n) in the maximum sensitivity         direction, which thus corresponds to the value of the         directivity diagram for the direction (θ_(n) ^(α),φ_(n)         ^(α))(f); and     -   parameters d_(n)(f) that are representative of the directivity         diagram of each sensor, which allows the directivity of each         sensor to be described by a model comprising a combination of         omnidirectional and bi-directional diagrams oriented in the         direction α_(n)(f), using the following directivity model:         1−d _(n)(f)+d _(n)(f)cos(α_(n)(f).(θ, φ))         wherein α_(n)(f).(θ, φ) designates the scalar product between         the directions α_(n)(f) and (θ, φ).

This parameter d_(n)(f) may be determined using standard methods for estimating parameters, for example by applying a method of least squares that provides the value d_(n)(f), which minimises the error between the real directivity diagram and the modelled directivity diagram.

Advantageously, the calibration step 20 also allows the parameter σ² _(n)(f) which corresponds to the power spectral density of the background noise of the sensors, to be determined. The signal issued by the sensor 2 _(n) is thus picked up during this step 20, in the absence of an acoustic field. The parameter σ² _(n)(f) is determined using methods for estimating power spectral density, such as the so-called periodogram method, for example.

Depending on the embodiments, all or some of sub-steps 22 to 28 are repeated, in order, for example, to allow a plurality of types of parameters to be determined, wherein some sub-steps may be common to the determination of various types of parameters.

The calibration step 20 may also be carried out using means other than those described, such as direct measuring means—for example, using means for optically measuring the position of each elementary sensor 2 _(n) relative to the centre 4 of the acquisition means 1.

Furthermore, the calibration step 20 may carry out a simulation, using a computer, for example, of signals that are representative of the acquisition capacities of the elementary sensors 2 _(n).

It would therefore appear that this calibration step 20 allows all or some of the parameters that are representative of the structural and/or electro-acoustic characteristics of the acquisition means 1, which are used during the step 30 involving the determination of the encoding filters, to be determined.

FIG. 5 illustrates a flow chart of an embodiment of the step 30 involving the determination of the encoding filters.

The step 30 comprises a sub-step 32 that involves the determination of a matrix B that is representative of the acquisition capacities of the acquisition means 1 or sampling matrix.

In the described embodiment, the matrix B is determined from the parameters {right arrow over (x)}_(n), H_(n)(f), d_(n)(f), α_(n)(f) and B_(n,l,m)(f) and is a matrix of size N over (L(f)+1)², comprising elements B_(n,l,m)(f), the index n designating the row n, and the indices (l,m) designating the column l²+l+m. The matrix B therefore has the following form: $\left\lbrack \left. \quad\begin{matrix} {B_{1,0,0}(f)} & {B_{1,1,{- 1}}(f)} & {B_{1,1,0}(f)} & {B_{1,1,1}(f)} & \ldots & {B_{1,L,{- L}}(f)} & \ldots & {B_{1,L,0}(f)} & \ldots & {B_{1,L,L}(f)} \\ {B_{2,0,0}(f)} & {B_{2,1,{- 1}}(f)} & {B_{2,1,0}(f)} & {B_{2,1,1}(f)} & \ldots & {B_{2,L,{- L}}(f)} & \ldots & {B_{2,L,0}(f)} & \ldots & {B_{2,L,L}(f)} \\ \vdots & \vdots & \vdots & \vdots & \quad & \vdots & \quad & \vdots & \quad & \vdots \\ {B_{N,0,0}(f)} & {B_{N,1,{- 1}}(f)} & {B_{N,1,0}(f)} & {B_{N,1,1}(f)} & \ldots & {B_{N,L,{- L}}(f)} & \ldots & {B_{N,L,0}(f)} & \ldots & {B_{N,L,L}(f)} \end{matrix} \right\rbrack \right.$

Specific elements of the matrix B may be determined directly during steps 10 or 20. The matrix B is then supplemented with elements determined from a modelling of the sensors.

In this embodiment, each sensor n is modelled by a point sensor placed in the position {right arrow over (x)}_(n), exhibiting a directivity composed of a combination of omnidirectional and bi-directional diagrams of proportion d_(n)(f), oriented in the direction α_(n)(f) and having a frequency response H_(n)(f).

The complementary elements B_(n,l,m)(f) are then determined according to the relationship: $\begin{matrix} {{B_{n,l,m}(f)} = {4\pi\quad{H_{n}(f)}j^{l} \times \left\{ {{\left( {1 - {d_{n}(f)}} \right){j_{l}\left( {kr}_{n} \right)}{y_{l}^{m}\left( {\theta_{n},\phi_{n}} \right)}} - {{{jd}_{n}(f)} \times}} \right.}} \\ {\left( {{{j_{l}^{*}\left( {kr}_{n} \right)}{y_{l}^{m}\left( {\theta_{n},\phi_{n}} \right)}u_{r}} - {\frac{j_{l}\left( {kr}_{n} \right)}{{kr}_{n}}{R_{l}^{m}\left( {\cos\quad\theta_{n}} \right)}{{trg}_{m}(\phi)}u_{\theta}} +} \right.} \\ {\left. \left. {\frac{{mj}_{l}\left( {kr}_{n} \right)}{{kr}_{n}\sin\quad\theta_{n}}{y_{l}^{- m}\left( {\theta_{n},\phi_{n}} \right)}u_{\theta}} \right) \right\}\quad{wherein}} \\ {{j_{l}^{*}\left( {kr}_{n} \right)} = \frac{{{lj}_{l - 1}\left( {kr}_{n} \right)} - {\left( {l + 1} \right){j_{l + 1}\left( {kr}_{n} \right)}}}{{2l} + 1}} \\ {{R_{l}^{m}\left( {\cos\quad\theta_{n}} \right)} = \left\{ \begin{matrix} {\sqrt{l\left( {l + 1} \right)}{P_{l}^{1}\left( {\cos\quad\theta_{n}} \right)}} & {{{when}\quad m} = 0} \\ {{\sqrt{\frac{\left( {l - m} \right)\left( {l + m + 1} \right)}{2}}{P_{l}^{m + 1}\left( {\cos\quad\theta_{n}} \right)}} -} & {{{when}\quad 1} \leq m \leq {l - 1}} \\ {{\sqrt{\frac{\left( {l + m} \right)\left( {l - m + 1} \right)}{2}}{P_{l}^{m - 1}\left( {\cos\quad\theta_{n}} \right)}} -} & \quad \\ {\sqrt{\frac{l}{2}}{P_{l}^{l - 1}\left( {\cos\quad\theta_{n}} \right)}} & {{{when}\quad m} = 1} \end{matrix} \right.} \end{matrix}$ and wherein u _(r)=sin θ_(n) sin θ_(n) ^(α)(f)cos(φ_(n)−φ_(n) ^(α)(f))+cos θ_(n) cos θ_(n) ^(α)(f) u _(θ)=cos θ_(n) sin θ_(n) ^(α)(f)cos(φ_(n)−φ_(n) ^(α)(f))−sin θ_(n) cos θ_(n) ^(α)(f) u _(φ)=sin θ_(n) ^(α)(f)sin(φ_(n) ^(α)(f)−φ_(n))

In the event of the sensors being oriented radially, the relationship admits a simpler expression: $\begin{matrix} {{B_{n,l,m}(f)} = {4\pi\quad{H_{n}(f)}j^{l}{y_{l}^{m}\left( {\theta_{n},\phi_{n}} \right)}}} \\ {\left( {{\left( {1 - {d_{n}(f)}} \right){j_{l}\left( {kr}_{n} \right)}} - {{{jd}_{n}(f)}\quad\frac{{{lj}_{l - 1}\left( {kr}_{n} \right)} - {\left( {l + 1} \right){j_{l + 1}\left( {kr}_{n} \right)}}}{{2l} + 1}}} \right)} \end{matrix}$

The step 30 then comprises a sub-step 34 involving the determination of an intercorrelation matrix A that is representative of the similarity between the signals c₁ to c_(N) issued by the sensors 2 ₁ to 2 _(N), owing to the fact that these sensors 2 ₁ to 2 _(N) carry out measurements on a single acoustic field P. The matrix A is determined from the sampling matrix B. A is a matrix of size N over N, obtained by means of the relationship: A=B B^(T)

Advantageously, the matrix A is determined more precisely using a matrix B that is supplemented up to an order L₂, according to the method of the preceding step.

Since the matrix A may be expressed solely as a function of the matrix B, the sub-step 34 involving the determination of the intercorrelation matrix A may be considered as an intermediate calculation step, and may thus be incorporated into another sub-step of the step 30.

The step 30 then comprises a sub-step 36 involving the determination of an encoding matrix E(f) that is representative of the encoding filters for a given frequency. The matrix E(f) is determined from the matrices A and B and from the parameters L(f), H(f), {(l_(k),m_(k))}(f) and σ_(n) ²(f). The matrix E(f) is a matrix of size (L(f)+1)² over N, comprising elements E_(l,m,n)(f), the indices (l,m) designating the row l²+l+m, and the index n designating the column n. The matrix E(f) therefore has the following form: $\left\lbrack \left. \quad\begin{matrix} {E_{0,0,1}(f)} & {E_{0,0,2}(f)} & \ldots & {E_{0,0,N}(f)} \\ {E_{1,{- 1},1}(f)} & {E_{1,{- 1},2}(f)} & \ldots & {E_{1,{- 1},N}(f)} \\ {E_{1,0,1}(f)} & {E_{1,0,2}(f)} & \ldots & {E_{1,0,N}(f)} \\ {E_{1,1,1}(f)} & {E_{1,1,2}(f)} & \ldots & {E_{1,1,N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {E_{L,{- L},1}(f)} & {E_{L,{- L},2}(f)} & \ldots & {E_{L,{- L},N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {E_{L,0,1}(f)} & {E_{L,0,2}(f)} & \ldots & {E_{L,0,N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {E_{L,L,1}(f)} & {E_{L,L,2}(f)} & \ldots & {E_{L,L,N}(f)} \end{matrix} \right\rbrack \right.$

The matrix E(f) is determined row by row. For each operating frequency f, each row E_(l,m) of index (l,m) of the matrix E(f) assumes the following form: [E_(l,m,1)(f)E_(l,m,2)(f) . . . E_(l,m,N)(f)]

The elements E_(l,m,n)(f) of the row E_(l,m) are obtained by the following expressions:

-   -   if (l,m) belongs to the list {(l_(k),m_(k))}(f), then:         E _(l,m)+μ(f)B _(l,m) ^(T)((μ(f)−λ)A+(1−μ(f)Σ_(N))⁻¹         wherein λ confirms the relationship:         (μ(f))² B _(l,m) ^(T)((μ(f)−λ)A+(1−μ(f))Σ_(N))⁻¹         A(μ(f)−λ)A+(1−μ(f))Σ_(N))⁻¹ B _(l,m)=1         and wherein the value of λ is determined using analytical or         numerical methods for investigating equation roots, optionally         using methods of matrix diagonalisation; and     -   if (l,m) does not belong to the list {(l_(k),m_(k))}(f), then:         E _(l,m)=μ(f)B _(l,m) ^(T)((μ(f)A+(1−μ(f))Σ_(N))⁻¹

In these expressions, B_(l,m) is the column (l,m) of the matrix B and Σ_(N) is a diagonal matrix of size N over N, which is representative of the background noise of the sensors, wherein the element n of the diagonal is σ_(n) ²(f).

The sub-steps 32, 34 and 36 involving the determination of the matrices A, B and E(f) are repeated for each operating frequency f.

Of course, in simplified embodiments, the parameters are frequency-independent, and the sub-steps 32, 34 and 36 are only carried out once. The sub-step 36 then allows directly the determination of a frequency-independent matrix E.

During a subsequent sub-step 38, parameters FD that are representative of the encoding filters are determined from the matrix E(f). Each element E_(l,m,n)(f) of the matrix E(f) represents the frequency response of an encoding filter. Each encoding filter may be described by the parameters FD, in different forms.

If, for example, the parameters FD that are representative of the filters E_(l,m,n)(f) are:

-   -   frequency responses, the parameters FD are then directly the         E_(l,m,n)(f) calculated for specific frequencies f;     -   finite impulse responses c_(l,m,n)(t) calculated by inverse         Fourier transformation of E_(l,m,n)(f), each impulse response         c_(l,m,n)(t) is sampled, then truncated to a suitable length for         each response; and     -   recursive filter coefficients with infinite impulse responses         calculated from E_(l,m,n)(f) using adaptation methods.

The step 30 involving the determination of the encoding filters thus issues parameters FD describing encoding filters that are representative of at least the structural and/or elctro-acoustic capacities of the acquisition means 1.

In particular, these filters are representative of the following characteristics:

-   -   position of the sensors 2 ₁ to 2 _(N);     -   intrinsic electro-acoustic characteristics of the sensors 2 ₁ to         2 _(N), in particular power spectral density of the background         noise and acquisition capacities of the acoustic field; and     -   optimisation strategies, in particular the compromise between         spatial faithfulness of acquisition of the acoustic field and         minimisation of the background noise produced by the sensors.

FIG. 6 illustrates in detail an embodiment of the step 40 involving the processing of the measurement signals issued by the acquisition means 1, by applying encoding filters to these signals and by adding the filtered signals.

In the step 40, the coefficients {circumflex over (p)}_(l,m)(t) that are representative of the acoustic field P are deduced from the signals c₁ to c_(N) derived from the elementary sensors 2 ₁ to 2 _(N), by applying the frequency-response encoding filters E_(l,m,n)(f) in the following manner: ${{\hat{P}}_{l,m}(f)} = {\sum\limits_{n = 1}^{N}{{E_{l,m,n}(f)}{C_{n}(f)}}}$ wherein {circumflex over (P)}_(l,m)(f) is the Fourier transform of {circumflex over (p)}_(l,m)(t) and C_(n)(f) is the Fourier transform of c_(n)(t).

The example described the case of filtering by finite impulse response. This filtering requires the determination, initially, of a parameter T_(n,l,m), corresponding to the suitable number of samples for each response e_(n,l,m)(t), which results in the following convolution expression: ${{\hat{P}}_{l,m}\lbrack t\rbrack} = {\sum\limits_{n = 1}^{N}{\sum\limits_{\tau = 0}^{T_{n,l,m} - 1}{{e_{n,l,m}\lbrack\tau\rbrack}{c_{n}\left\lbrack {t - \tau} \right\rbrack}}}}$

These coefficients {circumflex over (p)}_(l,m) are a finite number of coefficients that are representative over time and in the three-dimensional space of the acoustic field, and form a faithful representation of this acoustic field.

Depending on the nature of the parameters FD, other filtering processes by E_(l,m,n)(f) may be carried out according to various filtering methods, such as, for example:

-   -   if the parameters FD provide the frequency responses         E_(l,m,n)(f) directly, the filtering is carried out using         filtering methods in the frequency domain, such as block         convolution processes, for example;     -   if the parameters FD provide the finite impulse response         c_(l,m,n)(t), the filtering is carried out in the time domain by         convolution; and     -   if the parameters FD provide the coefficients of a recursive         filter with infinite impulse response, the filtering is carried         out in the time domain by means of the recurrence relation.

It would therefore appear that the invention allows an acoustic field to be represented faithfully, by means of a representation that is substantially independent of the characteristics of the acquisition means, in the form of Fourier-Bessel coefficients.

Moreover, as previously stated, the method of the invention may be carried out in simplified embodiments.

If, for example, all of the sensors 2 ₁ to 2 _(N) are substantially omnidirectional and substantially identical in terms of sensitivity and level of background noise, the method of the invention may be carried out solely on the basis of knowledge of the parameters {right arrow over (x)}_(n) that are representative of the position of the sensors 2 _(n) relative to the centre 4 of the acquisition means 1, and of the parameters μ and L, which relate to the optimisation strategy.

Moreover, in this simplified embodiment, the parameters are considered to be frequency-independent.

Using these parameters, the matrices A and B are thus calculated simultaneously or sequentially in any order during the sub-steps 32 and 34.

The elements B_(n,l,m)(f) of the matrix B are then organised in the following manner: $\left\lbrack \left. \quad\begin{matrix} {B_{1,0,0}(f)} & {B_{1,1,{- 1}}(f)} & {B_{1,1,0}(f)} & {B_{1,1,1}(f)} & \ldots & {B_{1,L,{- L}}(f)} & \ldots & {B_{1,L,0}(f)} & \ldots & {B_{1,L,L}(f)} \\ {B_{2,0,0}(f)} & {B_{2,1,{- 1}}(f)} & {B_{2,1,0}(f)} & {B_{2,1,1}(f)} & \ldots & {B_{2,L,{- L}}(f)} & \ldots & {B_{2,L,0}(f)} & \ldots & {B_{2,L,L}(f)} \\ \vdots & \vdots & \vdots & \vdots & \quad & \vdots & \quad & \vdots & \quad & \vdots \\ {B_{N,0,0}(f)} & {B_{N,1,{- 1}}(f)} & {B_{N,1,0}(f)} & {B_{N,1,1}(f)} & \ldots & {B_{N,L,{- L}}(f)} & \ldots & {B_{N,L,0}(f)} & \ldots & {B_{N,L,L}(f)} \end{matrix} \right\rbrack \right.$ wherein B _(n,l,m)(f)=4πj ^(l) j _(l)(kr _(n))y _(l) ^(m)(θ_(n),φ_(n))

Similarly, the elements A_(n1,n2)(f) of the matrix A are then organised in the following manner: $\left\lbrack \left. \quad\begin{matrix} {A_{1,1}(f)} & {A_{1,2}(f)} & \ldots & {A_{1,N}(f)} \\ {A_{2,1}(f)} & {A_{2,2}(f)} & \ldots & {A_{2,N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {A_{N,1}(f)} & {A_{N,2}(f)} & \ldots & {A_{N,N}(f)} \end{matrix} \right\rbrack \right.$

In this embodiment, the matrix A is obtained from the matrix B by means of the relationship: A=B B^(T)

Advantageously, the elements A_(n1,n2)(f) of the matrix A are determined with greater precision by means of the relationship: $\begin{matrix} {{A_{n_{1},n_{2}}(f)} = {4\pi{\sum\limits_{l = 0}^{L_{2}}{\left( {{2l} + 1} \right){j_{l}\left( {kr}_{n_{1}} \right)}{j_{l}\left( {kr}_{n_{2}} \right)}}}}} \\ {P_{l}\left( {{\cos\quad\theta_{n_{1}}\cos\quad\theta_{n_{2}}} + {\sin\quad\theta_{n_{1}}\sin\quad\theta_{n_{2}}{\cos\left( {\phi_{n_{1}} - \phi_{n_{2}}} \right)}}} \right)} \end{matrix}$ wherein L₂ is the order in which the determination of the matrix A is conducted and is an integer greater than L. The greater the value selected for L₂, the more precise, but longer, the calculation of the A_(n1,n2)(f) will be.

In the sub-step 36, the encoding matrix E which is representative of the encoding filters, is determined from the matrices A and B and the parameter μ according to the expression: E=μB ^(T)(μA+(1−μ)I _(N))⁻¹

The elements E_(l,m,n)(f) of the matrix E are organised in the following manner: $\begin{bmatrix} {\quad{E_{0,0,1}(f)}} & {E_{0,0,2}(f)} & \cdots & {E_{0,0,N}(f)} \\ {E_{1,{- 1},1}(f)} & {E_{1,{- 1},2}(f)} & \cdots & {E_{1,{- 1},N}(f)} \\ {\quad{E_{1,0,1}(f)}} & {E_{1,0,2}(f)} & \cdots & {E_{1,0,N}(f)} \\ {\quad{E_{1,1,1}(f)}} & {E_{1,1,2}(f)} & \cdots & {E_{1,1,N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {E_{L,{- L},1}(f)} & {E_{L,{- L},2}(f)} & \cdots & {E_{L,{- L},N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {\quad{E_{L,0,1}(f)}} & {E_{L,0,2}(f)} & \cdots & {E_{L,0,N}(f)} \\ \vdots & \vdots & \quad & \vdots \\ {\quad{E_{L,L,1}(f)}} & {E_{L,L,2}(f)} & \cdots & {E_{L,L,N}(f)} \end{bmatrix}\quad$

The sub-steps 32, 34 and 36 involving the determination of the matrices A and B, then E are repeated for all of the operating frequencies f

Each element E_(l,m,n)(f) corresponds to an encoding filter that incorporates the spatial distribution of the sensors 2 _(n) and also the optimisation strategy.

In the phase 40, the signals c₁ to c_(N) derived from the sensors 2 ₁ to 2 _(N) are filtered using encoding filters described by the parameters FD. Each coefficient {circumflex over (p)}_(l,m)(t) issued is deduced from signals c₁ to c_(N) by applying filters in the following manner: ${{\hat{P}}_{l,m}(f)} = {\sum\limits_{n = 1}^{N}{{E_{l,m,n}(f)}{C_{n}(f)}}}$ wherein {circumflex over (P)}_(l,m)(f) is the Fourier transform of {circumflex over (p)}_(l,m)(t), and C_(n)(f) is the Fourier transform of c_(n)(t).

In this embodiment, the coefficients {circumflex over (p)}_(l,m)(t) are determined using filtering methods in the frequency domain, such as block convolution methods, for example.

The representation of the acoustic field therefore takes into consideration the position of the sensors and the selected optimisation parameters and constitutes a faithful estimate of the acoustic field.

FIG. 7 is a block diagram of a device that is suitable for carrying out the method of the invention.

In this figure, a device 50 for representing the acoustic field P is connected to the acquisition means 1, as described with reference to FIG. 2.

The device 50, or encoding device, is also connected at the input to means 60 for determining parameters that are representative of the structural and/or electro-acoustic characteristics of the acquisition means 1.

These means 60 comprise, in particular, means 62 for inputting parameters and calibration means 64, which are suitable for carrying out steps 10 and 20, respectively, of the method of the invention, as described above.

The encoding device 50 receives, from means 60 for determining the parameters, a plurality of parameters that are representative of the characteristics of the acquisition means 1 that are distributed between a signal CL for defining the structural characteristics and a signal CP for the parameterisation of the structural and/or electro-acoustic characteristics.

The device also receives parameters relating to representation strategies in a signal OS for optimising representation.

In these signals, the parameters are distributed in the following manner:

-   -   in the definition signal CL:         -   parameters {right arrow over (x)}_(n) that are             representative of the position of the sensor 2 _(n);     -   in the parameterisation signal CP:         -   parameters H_(n)(f) that are representative of the frequency             response of the sensor 2 _(n);         -   parameters d_(n)(f) that are representative of the             directivity diagram of the sensor 2 _(n);         -   parameters α_(n)(f) that are representative of the             orientation of the sensor 2 _(n);         -   parameters σ² _(n)(f) that are representative of the power             spectral density of the background noise of the sensor 2             _(n); and         -   parameters B_(n,l,m)(f) that are representative of the             acquisition capacities of the sensor 2 _(n); and     -   in the optimisation signal OS:         -   a parameter μ(f) specifying the compromise between the             faithfulness of representation of the acoustic field and             minimisation of the background noise produced by the             sensors;         -   a parameter L(f) specifying the order in which the             representation is conducted; and         -   a parameter {(l_(k),m_(k))}(I) that is representative of the             list of the coefficients, the power of which must be equal             to the power of the corresponding coefficient in the             acoustic field to be represented P.

Advantageously, this device 50 comprises means 51 for formatting input signals that are suitable for issuing, from signals c₁ to c_(N), a corresponding formatted signal SI.

For example, the means 51 comprise analogue-digital converters, amplifiers or even filtering systems.

The device 50 further comprises means 52 for determining the encoding filters, which means comprise a module 55 for calculating the sampling matrix B and a module 56 for calculating the intercorrelation matrix A, both of which are connected to a module 57 for calculating the encoding matrix E(f).

This encoding matrix E(f) is used by a module 58 for determining encoding filters that issues a signal S_(FD), which contains the parameters FD that are representative of the encoding filters.

This signal S_(FD) is used by a processing module 59 that applies the encoding filters to the signal SI in order to issue a signal SI_(FB) which comprises the Fourier-Bessel coefficients that are representative of the acoustic field P.

Optionally, the device 50 comprises a non-volatile memory in which the parameters that form the signal S_(FD), which have been determined previously, are stored.

For example, the acquisition means 1 are tested and calibrated by their manufacturer in order to provide directly a memory comprising all of the parameters of the signal S_(FD) that are to be incorporated into an encoding device in order to acquire the acoustic field P and to issue a faithful representation thereof.

Similarly, in a variant, this memory comprises only the matrices B and optionally A, and the device 50 comprises means for inputting the parameters forming the optimisation signal OS, in order to carry out the determination of the encoding matrix E(f) and the determination of the parameters FD that are representative of the encoding filters.

Other distributions between the various modules described may, of course, be envisaged, as required. 

1. Method for representing an acoustic field comprising a step involving the acquisition of measurement signals issued by acquisition means comprising one or more elementary sensors that are exposed to said acoustic field, wherein said method comprising: a step involving the determination of encoding filters that are representative of at least the structural characteristics of said acquisition means; and a step involving the processing of said measurement signals by applying said encoding filters to these signals in order to determine a finite number of coefficients representative over time and in the three-dimensional space of said acoustic field, said coefficients allowing a representation of said acoustic field to be obtained that is substantially independent of the characteristics of said acquisition means.
 2. Method according to claim 1, wherein said structural characteristics comprise at least position characteristics of said elementary sensors relative to a predetermined reference point of said acquisition means.
 3. Method according to claim 1, wherein said encoding filters are also representative of electro-acoustic characteristics of said acquisition means.
 4. Method according to claim 3, characterised in that said electro-acoustic characteristics comprise at least characteristics related to the intrinsic electro-acoustic acquisition capacities of said elementary sensors.
 5. Method according to claim 1, wherein said coefficients allow a representation of the acoustic field to be obtained are what are known as Fourier-Bessel coefficients and/or linear combinations of Fourier-Bessel coefficients.
 6. Method according to claim 1, wherein said step involve the determination of the encoding filters comprises: a sub-step involving the determination of a sampling matrix that is representative of the acquisition capacities of said acquisition means; a sub-step involving the determination of an intercorrelation matrix that is representative of the similarity between said measurement signals issued by the elementary sensors forming said acquisition means; and a sub-step involving the determination of an encoding matrix from said sampling matrix, said intercorrelation matrix and a parameter that is representative of a desired compromise between faithfulness of representation of the acoustic field and minimisation of the background noise caused by the acquisition means, which matrix is representative of said encoding filters.
 7. Method according to claim 6, wherein said sub-steps involve the determination of the matrices are carried out for a finite number of operating frequencies.
 8. Method according to claim 6, wherein said step involve the determination of the sampling matrix is carried out, for each of said elementary sensors forming said acquisition means, from: parameters that are representative of the position of said sensor relative to the centre of said acquisition means; and/or a finite number of coefficients that are representative of the acquisition capacities of said sensor.
 9. Method according to claim 8, wherein said step involve the determination of the sampling matrix is also carried out from at least one of the following parameters: parameters that are representative of the frequency responses of all or some of the sensors; parameters that are representative of the directivity diagrams of all or some of the sensors; parameters that are representative of the orientations of all or some of the sensors, i.e. of their maximum sensitivity direction; parameters, that are representative of the power spectral densities of the background noise of all or some of the sensors; a parameter specifying the order in which the representation is conducted; a parameter that is representative of a list of coefficients, the power of which must be equal to the power of the corresponding coefficient in the acoustic field to be represented.
 10. Method according to claim 1 further comprising a calibration step allowing all or some of the parameters used in said step involving the determination of the encoding filters, to be issued.
 11. Method according to claim 10, wherein said calibration step comprises, for at least one of said elementary sensors forming said acquisition means: a sub-step involving the acquisition of signals that are representative of the acquisition capacities of said at least one sensor; and a sub-step involving the determination of parameters representative of electro-acoustic and/or structural characteristics of said at least one sensor.
 12. Method according to claim 11, wherein said calibration step further comprises: a sub-step involving the emission of a specific acoustic field toward said at least one sensor, said acquisition sub-step corresponding to the acquisition of the signals issued by this sensor when it is exposed to said specific acoustic field; and a sub-step involving the modelling of said specific acoustic field in a finite number of coefficients, in order to allow said sub-step involving the determination of parameters that are representative of electro-acoustic and/or structural characteristics of the sensor to be carried out.
 13. Method according to claim 10, wherein said calibration step comprises a sub-step involving the reception of a finite number of signals that are representative of the electro-acoustic and structural characteristics of said sensors forming said acquisition means, which signals are used directly during said sub-step involving the determination of the electro-acoustic and/or structural characteristics of said acquisition means.
 14. Method according to claim 1, wherein it comprises an input step, allowing all or some of the parameters used during said step involving the determination of the encoding filters, to be determined.
 15. Computer programme comprising programme code instructions for implementing the steps of the method according to claim 1, when said programme is executed on a computer.
 16. Movable support of the type comprising at least one operation processor and a non-volatile memory element, wherein said memory comprises a programme comprising code instructions for implementing the steps of the method according to claim 1, when said processor executes said programme.
 17. Device for representing an acoustic field that is connectable to acquisition means comprising one or more elementary sensors issuing measurement signals when they are exposed to said acoustic field, wherein it comprises a module for processing the measurement signals by applying encoding filters that are representative of at least the structural characteristics of said acquisition means to these measurement signals, in order to issue a signal that comprises a finite number of coefficients representative over time and in the three-dimensional space of said acoustic field, said coefficients allowing a representation of said acoustic field to be obtained that is substantially independent of the characteristics of said acquisition means.
 18. Device according to claim 17, wherein said encoding filters are also representative of electro-acoustic characteristics of said acquisition means.
 19. Device according claim 17, wherein it further comprises means for determining said encoding filters that are representative of structural and/or electro-acoustic characteristics of said acquisition means.
 20. Device according to claim 19, wherein said means for determining encoding filters receive at the input at least one of the following parameters: parameters that are representative of the positions, relative to centre of said acquisition means, of all or some of the sensors; a finite number of coefficients that are representative of the acquisition capacities of all or some of the sensors; parameters that are representative of the frequency responses of all or some of the sensors; parameters that are representative of the directivity patterns of all or some of the sensors; parameters that are representative of the orientations of all or some of the sensors, i.e. of their maximum sensitivity direction; parameters that are representative of the power spectral densities of the background noise of all or some of the sensors; a parameter that is representative of the desired compromise between faithfulness of representation of the acoustic field and minimisation of the background noise caused by the acquisition means; a parameter specifying the order in which the encoding is conducted; a parameter that is representative of a list of coefficients, the power of which must be equal to the power of the corresponding coefficient in the acoustic field to be represented.
 21. Device according to claim 20, wherein it is associated with means for determining all or some of the parameters received by said means for determining the encoding filters, said means comprising at least one of the following elements: means for inputting parameters; and/or calibration means.
 22. Device according to claim 1, wherein it is associated with means for formatting said measurement signals, in order to issue a corresponding formatted signal. 